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Abstract 



^ . It has been recently shown that smah subsystems of finite quantum systems 



generically equilibrate. We extend these results to infinite-dimensional Hilbert 



Tij- \ spaces of field theories and matrix models. We consider a quench setup, where 



initial states are chosen from a microcanonical ensemble of finite energy in free 
theory, and then evolve with an arbitrary non- perturb at ive Hamiltonian. Given 



o 

ff^ ■ a dynamical assumption on the expectation value of particle number density, we 

prove that small subsystems reach equilibrium at the level of quantum wave- 
function, and with respect to all observables. The picture that emerges is that 

^ \ at higher energies, larger subsystems can reach equilibrium. For bosonic fields 

j^ I on a lattice, in the limit of large number of bosons per site, all subsystem smaller 

than half equilibrate. In the Hermitian matrix model, by contrast, this occurs 
in the limit of large energy per matrix element, emphasizing the importance 
of the 0{N'^) energy scale for the fast scrambling conjecture. Applying our 
techniques to continuum field theories on compact spaces, we show that the 
density matrix of small momentum-space observables equilibrate. Finally, we 
discuss the connection with scrambling, and provide a sufficient condition for a 
time-independent Hamiltonian to be a scrambler in terms of the entanglement 
entropy of its energy eigenstates. 
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1 Introduction 

It is generally accepted that most systems equilibrate even if they are perturbed far 
away from equilibrium. This is to say that the information about the initial state 



spreads out such that after some time all observables restricted to small subsystems 
become almost independent of time and the initial state. Our current understanding 
of equilibration in field theories remains at the level of perturbation theory, either at 
small coupling or large A^ expansion. Perturbation theory enables us to keep track 
of the time dependence of the expectation value of observables with a small number 
of field operators, what we refer to as small observables. However, the formalism 
quickly becomes cumbersome as we look at large observables with many field insertions. 
Knowledge of all observables is required to specify the quantum state of field theory. 
It is not clear whether the effective state one infers from the expectation value of small 
observables is a good description of the system. 

In this work, we address equilibration at the level of subsystems' density matrices 
non-perturbatively. We argue, quite generally, that the reduced density matrix of 
small subsystems becomes indistinguishable from the equilibrium state with respect 
to all observables living in the subsystem Hilbert space. The techniques we employ 
in this work apply to Hamiltonians with discrete spectrum as in lattice field theory, 
continuum field theories on compact manifolds and matrix models. Our results suggest 
that at higher energies larger subsystems can equilibrate and in the limit of infinite 
energy, on a lattice and in matrix models, any subsystem smaller than half the size of 
the system equilibrates. 

Our work is a generalization of quantum information techniques [H 121 [3] that have 
been recently used to discuss equilibration of finite systems to the case of infinite 
dimensional Hilbert spaces]^ Motivated by the fast scrambling conjecture [5], this 
generalization enables us to address the question of the equilibration of large subsystems 
in field theories and matrix models. We consider a quench setting where the initial 
state is chosen from a free theory, and evolved with an arbitrary non-perturbative 
interacting Hamiltonian. With some reasonable physical assumptions we prove that 
small subsystems always equilibrate. Furthermore, on a lattice of A^ sites, initial states 
of energy larger than 0{N) equilibrate on all subsystems larger than half. In the 
Hermitian matrix model the large-scale equilibration occurs for energies larger than 
0{N^). In continuum field theories, our results imply that there are scales /ii and 
/i2, where the infra-red and ultra-violet density matrices constructed by tracing out 
all momentum modes, respectively, higher than /xi and lower than fi2, are almost non- 
dynamical and time-independent. In a quench setting from an eigenspace of large 
energy A in a (i-dimensional free theory, we find that both /ii and /i2 are 0(A^/'^). 

In spite of the technical details of proofs, the equilibration results are based on a 
rather simple picture. The initial state of the system lives in a large Hilbert space and 
contains a significant amount of information. If the Hilbert space that the subsystem 



^For a discussion of equilibration in integrable systems see [3]. 



explores is small, then it can not contain all the information about the initial state 
and quickly loses track of it. An interesting question is what happens if one tries 
to encode a message in the initial state. We say that the message is scrambled if 
it cannot be retrieved from the state of any subsystem smaller than half the size of 
the system [3 E]. With this definition, scrambling is a strong form of equilibration 
in which the state of large subsystems become independent of initial states including 
atypical ones such as a tensor product states. Hamiltonians, which are very powerful 
in generating entanglement, leave almost no room in the Hilbert space of subsystems, 
even the large ones. We will discuss the implication of this for scrambling of time- 
independent Hamiltonians. 

In section 2, we start by defining what we mean by equilibration. Our method 
closely follows the arguments and approach to equilibration developped in the quan- 
tum information community over the past few years. Therefore, we devote section 

3 to a quick review of results in finite systems. Our main results appear in section 

4 which generalizes equilibration theorems to infinite dimensional Hilbert spaces in 
quench settings and includes as an example the equilibration in the following three 
systems: bosons on a lattice, the Hermitian matrix model and field theory in momen- 
tum space. All results up to section 5 hold for generic interactions. The dynamical 
input from the Hamiltonian appears in section 5 where we expand on the connection 
with scrambling. 



2 Basic setup 

The many-particle quantum systems we are interested in are described by a set of 
physical states living in a tensor product of Fock spaces: ip & J^ = J^i ® J^2 ® ■ ■ ■ J^n- 
The tower of states in Fock spaces J-^ are generated by the creation operators a]. In 
particular we consider three systems described in this manner: 

1. Bosons on a lattice where a] creates a boson at site i. 

2. An N X N Hermitian matrix model which is equivalent A^ interacting fermions. 
The creation operator a] increases the energy of i*^ fermion by one unit. 

3. Quantum fields in momentum space on a compact space where a|. creates a particle 
with momentum k. 

Note that for bosons on a lattice and quantum fields in momentum space, any ip E T 
is physical, whereas in the matrix model the physical states are only the ones invariant 
under gauge transformations. In the Hermitian matrix model we consider here, due to 



the fermionic nature of degrees of freedom, the physical states are ant i- symmetrized 
over A^ sites [7]. 

In physical problems of interest, there are typically global constraints on the Hilbert 
space of initial states. They are chosen from a subspace I-Lr of J-". Note that the 
restriction is imposed only on initial states and the wave-function can leave I-Lr as it 
evolves. We denote the subspace the time-evolved initial states explore by T-Lt- It is 
the subspace spanned by energy eigenstates with nonzero projection into I-Lr. 

Consider an initial physical state ip G Hr d J^. K^ ip evolves in time, if the system 
equilibrates with respect to subsystem S, the reduced density matrix on S becomes 
close to an equilibrium state p^g^^ . Any comparison of distances between physical 
states requires a choice of metric in the Hilbert space. Trace norm provides a metric 
with a natural operational meaning. If two states ps and p^^™ are close in trace 
norm, i.e. Hp^ — p^^"* ||i < e, then any measurement to distinguish them succeeds with 
probability at most ^ + eo For some choice of Hamiltonians and initial states it is 
possible that that ps{t) approaches its equilibrium state monotically in time [HI ITO] . 
We say ps equilibrates in a strong sense, if for any e > there exists a r such that for 
all t > T, ||ps(t) — p^^"* 111 < e. More generically, systems tend to equilibrate in a weak 
sense: ps{t) becomes close to its equilibrium value, and spends most of its time in a 
small neighbourhood around it. We say ps equilibrates in a weak sense, if there exists 
an e <C 1 such that: 

(iiP5(t)-pr"iii)*<e, (2-1) 

where {cr{t))t = limr^oo j^ Jq <j{t)dt] see figure 1. 

Conceptually, the weak equilibration matches our intuition of equilibration better, 
since it leaves the possibility for ps to fluctuate away from equilibrium in a time- 
independent way. The bound in (12. ip is the statement that large fluctuations are rare. 

We would like to stay as generic as possible; therefore we focus on weak equilibra- 
tion. Following [1] , our strategy is to show first that the distance between a subsystem's 
density matrix and its time-averages u = \mvT^oo{ps{^))t is bounded by a constant: 

{\\ps{t)-uJs\\i)t<Tls. (2.2) 

Then, we show that for most initial states in the initial ensemble the constant rjs is 
small. This establishes the time-independence of the subsystem state. In principle us is 
a function of the initial state if). Next, we need to show that uJs{ip) is well-approximated 
by its average over all ip G I-Lr: 

{\\ujs{'iP)-{uJs)M'^<^'- (2-3) 



^For a basic reference on trace norm, and other commonly used measures of distance in Hilbert 
space see [5]. 




Figure 1: This figure schematically shows the difference between weak and strong 
equilibrations. As the subsystem density matrix evolves it comes close to its equilibrium 
value with respect to the trace norm. In weak equilibration (left-hand side) ps{t) spends 
most of its time in an e neighborhood of the equilibrium state which is independent of 
the initial state po, (pQ or xo- This is in contrast with strong equilibration (right-hand 
side) where equilibrium is approached monotically in time. 



This is the statement that the equilibrium state keeps almost no memory of the 
initial perturbation. However, this does not necessarily imply that the equilibrium state 
is well approximated by the Gibbs state. In cases where the asymptotic state is the 
Gibbs state we say that the system has thermalized. It is known that thermalization 
is not as ubiquitous in nature as equilibration. There are many known examples of 
systems and initial states that equilibrate but never thermalize [101 IH] . It is a trivial 
extension of the above definitions to introduce notions of weak or strong scrambling or 
thermalization. 

Another noteworthy comment is the difference between the equilibration of sub- 
systems as defined above and equilibration with respect to coarse-grained observables 
[21 [3l [12]. Coarse-grained observables need not be the result of only local measure- 
ments. In appendix [B] we briefly summarize the current understanding of equilibration 
for coarse-grained observables. 



3 Intuition from finite systems 



Recent developements in quantum information theory have made it possible to reformu- 
late the foundations of thermodynamics in a more precise way. The authors in [H [21 [3] , 
among many others, revisited the question of equilibration in the quantum world and 
proved the following three statements about weak equilibration in finite-dimensional 
systems that we refer to as fluctuation, typicality and universality theorems: 



Theorem 3.1 Fluctuations: Consider an arbitrary initial state ipo in a finite dimen- 
sional Hilbert space Hr C "H evolving with a Hamiltonian with non-degenerate energy 
gapsU Let u = {'ip{t))t be the time-averaged density matrix of the whole system, and 
deffioj) = ^j.(Ij2) be a measure of the effective dimension of the Hilbert space it explores. 
For any subsystem S we define: 



Vsii^o 
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deffiu) 



^eff 



where ds is the dimension of ps- Then, 

{\\psit,iJo) -cus{M\)t<Vs{^o)- (3.1) 

Therefore, for all initial states ipo, psit, ipo) approaches its time average in a weak sence 

ifvsM^:'^- 

Hereafter, we refer to r]s as the equilibration parameter. Roughly speaking, rjs 
bounds the size of fiucutations around equilibrium, we expect r]s to be small whenever 
the subsystem S is much smaller than the system. One way to make this intuition 
more percise is to ask: what is the probability of finding a large deff{u) if we pick a 
random ipo from H/?? The answer to this question is given by the typicality theorem 
which makes use of Levy's lemma |j: 

Theorem 3.2 Typicality: For a state ipQ chosen at random with uniform measure 
in Hr: 

Prob (d,ff{uj) < — j < 2e-'=^^, (3.2) 

where c = ^^ ^2x 10^^ 

Putting the first two theorems together we conclude that as long as the dimension 
of the initial ensemble's Hilbert space Hr is much larger than the dimension of a 
subsystem's Hilbert space, any interacting Hamiltonian evolves the subsystem towards 
equilibration. Naturally, we expects the equilibrium state to be almost independent 
of the initial state ipo- The universality theorem quantifies this by showing that the 
equilibrium states corresponding to all states in T-Cr are almost indistinguishable: 



•^For any equilibration to happen in Hilbert space it is essential that the degrees of freedom interact. 
One way to assure this is to restrict the problem to Hamiltonians with no energy gap Ei — Ej that is 
hugely degenerate. In a free theory E ~ Es + Eg, so for any arbitrary Eg and energy gap SEgin S, 
there is a gap SE in E. Note that even an arbitrarily small generic interaction lifts all degeneracies. 
The generalization of this result to the case of degenerate interacting Hamiltonians is discussed in 

m- 

*For completeness, we have included Levy's lemma in appendix \X\ 
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ps{t,ipo) ► usiipo) = {ps{t,i'o))t *^ Us = {uJs{i'o))i,o ► e~^^ 

time-independence initial-state independence Gibbs state 

Equilibration thermalization 

Figure 2: Information-theoretic approach to equihbration: First, the fluctuation and 
typicahty theorems estabhsh the time-independence of the subsystem density matrix 
by showing that ps(t) is close to its time-average for a typical state. Then, the time- 
averaged state is shown to be close to its initial ensemble average. 



Theorem 3.3 Universality: For a random state tpQ chosen with uniform measure in 

Prob(\\{u{ij,),{co)^,\\ > ^yj + ^j <2e-^'^''"' (3-3) 

where c' = ■^. Here (w)^o ^^ ^^^ state a;(^o) averaged over "H^jlj 



Figure 2 summarizes the approach discussed above. 



4 Equilibration in infinite-dimensional Hilbert spaces 

In systems with bosonic degrees of freedom, ps lives in an infinite dimensional Hilbert 
space. Equilibration in trace norm is the condition that ps at any time is indistin- 
guishable from its equilibrium value with respect to an infinite number of independent 
measurements. It is not apriori clear that such a strong condition holds in equilibrated 
field theories. This problem is manifest in the appearance of ds in the equilibration 
parameter. Subsystems can have infinite dimensions, and since the denominator in 
rjs only depends on the initial state and is finite when d^ is finite, the bound seems 
useless. 

What saves us is the intuition that in any physical equilibration process the support 
of the wave function on arbitrarily large energy eigenstates is small. Therefore, one 
should be able to truncate the Hilbert space by reasonable dynamical assumptions 
without missing out the relevant physics. As long as the errors caused by truncation 



^The tighter version of this bound introduced in [T] is crucial for our discussion of scrambling in 
section 5. However, this form suffices for our purposes here. 



are small, the finite dimensional calculations are reliable. When the equilibration 
parameter is small but truncation errors are not, it is only the set of observables which 
act on the truncated subspace that equilibrate. Significant deviations from equilibrium 
can be detected using operators that act beyond the truncated space. 

Cutting off the full Hilbert space in the interacting Hamiltonian basis at energy A 
does not guarantee that the reduced density matrix on S is finite-dimensional. However, 
if we instead truncate in the free Hamiltonian basis, not only does it regulate ds but 
also it allows us to estimate the equilibration parameter by counting partitions. 

In order to make the discussion more concrete, we choose to study the following 
"quench" problem jj Consider a free theory on J-" = ^-"1 ® ... ® J-'n with Hamiltonian 
Hfree = Si=i f^iHi- The Hilbert space is a direct sum "H = (Be'He, where TIe is the 
eigenspace corresponding to energy E spanned by e^; = {|ni, . . . , un) '■ E = J2i l^i'^i}- 
We choose our initial states from a subspace corresponding to energy A: ?/;o G I-Lr = 
HaU Initially dji is the cardinality of ca, and ds is the number of distinct \nkj , . . . , UkJ 
one finds in ca, by restricting to subsystem's degrees of freedom: S = {ki, . . . ,ks}. 

At t = we quench the system by turning on an interaction Hamiltonian, and study 
the relaxation of ps{t)- If we assume that [Hint, Hfree] = 0, it is clear that ip{t) never 
leaves "Ha ,and the problem remains finite-dimensional. Hence, the results of previous 
section apply: subsystem S equihbrates if (i| <C d^. Generically [Hint, Hfree] 7^ 0, and 
both ip and ps have infinite-dimensional Hilbert spaces available to them to explore. 
Since d^ only depends on the initial ensemble it remains finite, whereas ds is unbounded 
above and rjs diverges. This is an artifact of the infinite dimensionality of the Hilbert 
space. Let us define the projector to the finite-dimensional Hilbert space (BekaHe'- 

A 

^^IV^)= Yl \ni,...,nN){ni,...,nN\4j). (4.1) 

In a similar fashion, P^ denotes the projector to the initial ensemble "Ha- We expect 
that there exists a truncation such that Psit) = tr s{P^^{t)P^) approximates psit) 
well. Here 5" denotes the complement of S. The triangle inequality tells us: 

||p^(t) -C.5II < ||pf(t) -cfll + \\{ps{t) - cos) - {pl{t)-col)\\ (4.2) 

We call the second term on the right hand side the truncation error, and denote it by 
£^. The trick is to choose a cut-off A such that both the equilibration parameter in 
the truncated subspace rjg, and the errors caused by truncation are small. For this 



^Typically the term quench is used when the initial state is the ground state of an initial Hamil- 
tonian. Here, the initial state is chosen from some large energy sector of the free theory. 

^Choosing ipQ £ ®e=a-s^e with some small (S/A will not change any of our general conclusions. 
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procedure to work, we need to make a dynamical assumption that limits the support 
of wave-function on I-Le for E > h., and quantifies the truncation errors. Here, we 
employ a generalization of the method in [H] that uses the mean photon number in 
quantum optics to quantify errors caused by truncation to a finite-dimensional Hilbert 
space. The following theorem provides upper bounds on both terms on the right hand 
side of g2D: 

Theorem 4.1 Equilibration in infinite-dimensional "H: Consider an arbitrary 
state ipo, chosen at random from a subspace corresponding to energy A in the Hilbert 
spcae of a free theory. For any subsystem S, and projectors P^ and P^ defined as in 
([^.i|), we define: 



vUi^o) - ' ^^^^ 



deffiu) 

where dg is the dimension of the projected Hilbert space available to S. After a quench 
by an arbitrary interaction Hamiltonian: 

(i) Fluctuations: 

{\\ps{t)''-ujl\\),<4. (4.3) 

(ii) Truncation error: 



£^ ^ {\\{ps{t)-ujs) - (Pt(t) -u;^)||). < Q.Mt^EllA (4.4) 



A 



(Hi) Typicality: 



Prob [d^fM < ^) < 2e-^v^, (4.5) 



where c = &- ^ 2 x 10"^ 

iZTV 

(iv) Universality: 



Prob I \\u;|{^Po) - {u^lUJ > U^ + H < 2e-'='"'^^ (4.6) 

S^ = J2(^\^\^)tT[tTsiP^\k){k\P^n (4.7) 

k 

where d = ^. Here {uj)^^^ is the state u^^ipo) averaged over l-i/^, and \k) are 
energy eigenstate of the interacting Hamiltonian. 



The above theorem iinphes that the state of subsystem S equihbrates for almost all 
initial states ipo G Ha, if there exist a cut-off A such that 77^ ^ 1 and tr (uHfree) ^ A. 
The proof appears in appendix O 

The dynamical condition to meet to justify a truncation with P^ is that the time- 
averaged expectation value of the free part of the Hamiltonian, which was initially A, 
is much smaller than the cut-off scale Ajf| For systems with A^ degrees of freedom, it 
is more convenient to define new variables p, m, rh and A: 

m = A/NP 

m = tr {LuHfree)/NP 

A = A/A, (4.8) 

that remain finite in the thermodynamic limit A^ — )■ 00. In terms of these new variables, 
the truncation error S^ < 6w ;^. We devote the rest of this section to discussing the 
implications of theorem 14.11 in three different physical systems. 

4.1 Bosons on a lattice 

As our first example of equilibration in multi-particle quantum systems we consider 
bosons on a lattice with A^ sites. The Hamiltonian of the free theory, in appropriate 
units, is simply the total boson number: 

H = fiN = n^ alai, (4.9) 

i 

where a, and a] are the creation and annihilation operators defined at each site on the 
lattice. Setting /i = 1, the energy spectrum is given by all non-negative integers A, 
each of which appear with degeneracy d^ = { ^_^ ) , which is the number of weak 
compositions of A in A^ partsjj We restict the ensemble of initial states to Ha- The 
effective dimension of ipo only depends on the initial ensemble, and for almost all ipo 
is larger than cIa/^- Whereas, the dimension of the subspace available to subsystem 
S after introducing the cut-off A is the number of distinct compositions one finds by 
restrciting to a sequence of length s = IS"! in weak compositions of A in A^ parts: 

We are interested in investigating how subsystems of different sizes equilibrate in the 
thermodynamic limit A^ — )■ 00: 



^Roughly speaking, this is a constraint on the number density of particles. 

^A weak composition of A in N parts is a sequence of N non-negative integers which sum to A. 
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Finite-site subystems: The density matrix of subsystems of size s = 0(1) con- 
tains information about small observables of the sort {(j)"^^{xi) . . .(j)'^''{xs)). Let us 
assume for the moment that the number of bosons per site is initially finite and re- 
mains so at all times, i.e. p < 1. From f l4.4p we find that the error caused by truncation 
by A = iV^ is negligible for any e > 0. Plugging these into the expression for the equi- 
libration parameter gives 



A ^ N'A i{mA + s-iy?{N-iy.A\ 



''^^^rV (iV^A)!2(A + iV-l)! ^^-^^^ 



Using the Stirling approximatioiiLj we find that: 

^ f f{s) Ns{p+e)+p/A-{l~p)mNP/2 if Q < p < 1 

for some 0(1) functions f{s) and g{s). Hence, in the thermodynamic limit A^ ^ 1 the 
density matrix of all finite subsystems equilibrate. 

A small but finite fraction of all sites: Consider subsystems of size s/N order 
one but small, with a finite initial boson number density per site: m = 0(1). Applying 
the Stirling approximation one finds that the equilibration parameter is bounded by 



V s J (1 + m) I \l + ml \s'^ 



„^<.if!^^!^V'^-,^] f^y"""7^ 




^ '™ivw!y-'- . (4,13) 



for some 0(1) constant c. For any initial energy per site m and small error S^, there 
exists an sq found from 

( ^^J =(l + m) (4.14) 

such that any subsystem smaller than sq equilibrates. 

Large subsystems: A careful look at equation (14.141) shows that for any fixed 
A, So{m) is an increasing function of m with the asymptotic value limm-j.oo Sq = A/2; 
see figure 14.11 . It is a curious fact that for large enough initial boson number density, 
not only small subsystems but also all subsystems smaller than half the size of the 
system equilibrate. Intuitively, this is due to the exponential growth in density of 



lOi <^ n\ ^ _e_ 

~ V2irn(ri/e)'» ~ 27r ' 
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Figure 3: The size of the largest subsystem that provably equihbrates plotted as a 
function of the initial average boson number per site. The error parameter is taken to 
be £^''* < 6 X \{^~^ -Jmjvfi. In the limit of infinite initial energy any subsystem smaller 
than half the size of the whole system equilibrates. 



states as a function of energy which makes the fraction of the Hilbert space available 
to S vanishingly small for A^ ^ 1. In fact, for any A > 0{N), in the limit of A^ — )■ oo, 
both the error and the equilibration parameter vanish for any subsystem smaller than 
half. 



4.2 Hermitian matrix model 

For bosons on a lattice we found that states with energies A > 0{N) equilibrate on 
all subsystems smaller than half. However, we do not expect such a sector to survive 
in continuum field theories. An example where one is interested in studying such high 
energy states is the matrix model. Consider the Hermitian one-matrix model with a 
free Hamiltonian 



(- 



H = tT iX^ + m^X' 



(4.15) 



where X is an A^ x A^ matrix. The Hamiltonian is U{N) invariant, a freedom which 
can be partially fixed by diagnalizing the matrix. In the diagnal form, each eigenvalue 
Aj becomes a simple harmonic oscillator. The change of variables in the quantum 
partition function gives a Van der Mond determinant which can be absorbed by a 
redefinition of eigenvalues [71 |T5] . In terms of redefined eigenvalues, the wave-function 
ip{Xi, . . . , Xn) is a tensor product of A^ harmonic oscillators states, anti-symmetrized 
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under the permutation group of eigenvalues: 



ip{ni, 



.riN 



W! 



det 



/|ni(Ai)) |ni(A2)) 

'^2(Al)) \n2{\2)) 



|^l(AAr))\ 



(4.16) 



\|?^Af(Ai)) |r2Ar(A2)) •■■ \nN{XN))/ 



where \ni{Xj)) is the n* excited state of the j eigenvalue and multiplication inside 
the determinant is tensor product of states. The fermionic nature of the degrees of 
freedom implies that Ui is ordered as un > riN^i > ... > rii. The ground state has 
energy 

1\ N^ 



^•°'-|(*4) 



UJ 



-00. 



(4.17) 



We will set a; = 1 and measure energy in units of u. The excitations are described by 
an ordered set of positive numbers r^v > tat-i > • • • ?"i > where each r^ corresponds 
to the number of times the i^^ energy level is excited over its ground state, i.e. 



r,; 



E,, - E. 



(0) 



(4.18) 



If A = £■ — E'^'^\ then to each total energy E corresponds a sub-Hilbert space "Ha 
spanned by degenerate eigenstates of that energy level. The dimension of "Ha is 
P{A,N), the number of paritions of A into at most A^ parts. To each of these states 
one can associate a Young diagrams with rjv_j+i boxes in the i^'^ row from the top 
[T5| [T6] . Then, P{A,N) counts the number of Young diagrams with A boxes and 
maximum A^ rows. 

Consider a subsystem of s eigenvalues. The dimension of dg is bounded above by 
the number of distinct subdiagrams one finds by picking s rows from Young diagrams 
corresponding to P(A, A^), and discarding the rest of the rows. Sub-diagrams have any 
number of boxes from one to A, but have at most s rows. Therefore, 



4<J2PiE,s)<AP{A,s). 



E=l 



Then, the equilibration parameter satisfies 



^s<Ai 



P(A,s)2 

p{A,Ny 



(4.19) 



(4.20) 



By studying the large A^ asymptotics of restricted partitions functions P{N, s) we 
finfl 

^^See appendix [Pj 
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A finite set of eigenvalues: For s = 0(1) and p > 0, the equilibration parameter 
is small: 

r^l < cA«e-^/2(^^'''""'"'') (4.21) 

Choosing A to be any polynomial of finite order larger than one in A keeps both the 
truncation error and the equilibration parameter small. Hence, in the thermodynamic 
limit, small subsystems equilibrate. 

Large subsystems: One particularly interesting regime is A = 0{N'^). In this 
regime, we find that for large A^ and s/N = 0(1), the equilibration parameter is 
bounded by 

r/^<cyiVe^("''"■^)^/^ (4.22) 

for an order one constant c. The exponent on the right hand side is given by: 



(/3c^ - ^ log(l - e-'^o) - (/3 - 2mlog(l - e'^) 



f{s,m,X) =2 (/?,-- 
where /3 and /3c are found by solving 



(4.23) 



m 



iog(i - e-^y 



s^ log(l — e~^=)^ 



(4.24) 



For fixed m and arbitrary small error S^ < Q^Jrh/^mX) all subsystems of size smaller 
than sqN equilibrate, where sq is the solution to /(sq, m, A) = 0. In the limit of a large 
initial boson number m ^ 1, and small errors £^'*' ^ 1, it is possible to simplify the 
expressions in (14.231) to obtain: 

m ~ ( -2^ j (4.25) 

The interesting fact is that by increasing m it becomes possible for larger and larger 
subsystems to equilibrate. In the energy range A > 0{N'^) we find that all subsystems 
smaller than half the size of the system equilibrate with a vanishing error. Note that in 
the Hermitian matrix model, the regime relevant for equilibration of large subsystems 
has an initial energy per degree of freedom that scales with A^, as opposed to bosons 
on a lattice where it was finite. 
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Figure 4: The size of the largest subsystem that provably equihbrates plotted as a 
function of the average initial energy per matrix element; A/N"^. The error parameter 
is suppressed hy E"^ < 6 x lQ~^^Jm/m. In the limit of infinite initial energy any 
subsystem smaller than half the size of the whole system equilibrates. 



4.3 Field theories in momentum space 

Our discussion so far has been restricted to discrete models. In continuum field theories 
the Hilbert space of free theory has a natural decomposition in terms of momentum 
modes. On a compact space and in the absence of any internal symmetries, the Hilbert 
space is a tensor product of infinite Fock spaces J^k corresponding to momentum modes 
k. There is no analogue of parameter N in momentum space, so we only consider high 
energy ensembles, i.e. A ^ 1. 

1+1-dimensional field theories 

For simplicity, let us start with a scalar field in a one-dimenstional box and normalize 
such that /c's are positive integerso The Hamiltonian of a free massless theory has the 
form 



H = y^^kalak- 



(4.26) 



The dimension of the energy eigenspace A is given by the number of partitions of A, 
which according to the Hardy- Ramanuj an expansion scales as 



^A = P{A) ^ e"V2K75_ 



(4.27) 



12 



The generalization to fields with non-zero spin is straightforward. 
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On the other hand, the dimension of the Hilbert space available for a set of momentum 
modes S = {ki, . . . , ks} with a momentum cut-off A satisfies 



d^<J2DiE;S), (4.28) 



E=l 

and D{E; S) is the denumerant; the number of ways one can partition the number E 
as aiki + . . . + Qskg. 

Finite momentum modes: The large E asymptotic of denumerant for finite s is 
given bjo 

It is clear from f l4.29p that for large A, dg is bounded by a polynomial of degree s in 
A. As a result, for any A polynomial in A of a finite order larger than one, both the 
equilibration parameter and the truncation error are small. The theorem 14.11 implies 
that ps equilibrates for any finite set of momentum modes. In other words, the results of 
all field theory measurements and observables that involve a finite number of k modes, 
almost at all times after the quench, are well approximated by those of a stationary 
state, conditioned that whatever the non-perturbative dynamics is, the expectation 
value of kinetic energy does not grow indefinitely in time. 

Large subsystems: An interesting observation is that for any subsystem size s, dg 
is at its maximum value when S is the "infra-red" density matrix; i.e. S = {k : \k\ < 
g||i4| rjj^g reason is that for all S and E the denumerant satisfies D{E; S) < P{E, s), 
where P{E, s) is the number of partitions of E into at most s parts with the equality 
holding when S = {k : \k\ < s}. l£| 

Let pjj^ and p^^ denote, respectively, the "infra-red" and "ultra-violet" density 
matrices corresponding to the set of modes S = {k : k < p} and S = {k : k > p}. For 
a finite energy state, the expecation is that we do not lose details of the dynamics by 
tracing out large momentum modes. In this context of equilibration, we ask what are 
the largest and smallest) p such that, respectively, pj/^ and p^^ are time-independent, 
and contains no information about the dynamics or the initial state. 



^^See appendix [Pl for details. 

^ Note that in a theory with non-perturbative interactions, fow-momentum modes can have overlap 
with large energy eigenstates. Therefore, the reduced density matrix ps is not necessarily capturing the 
low-energy dynamics. However, the constraint on the time-average of the free Hamiltonian expectation 
value implies that this overlap is small for eigenstates of energy larger than A. 

^^We are using the following well-known result in number theory: the number of partitions of E 
into at most m parts is equal to the number of partitions of E with each part at most m. 
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The dimension of the Hilbert space available to p^J^ is 

A 

<^<E^(^'^)- (4.30) 



E=l 



The asymptotic behaviour of P{E, /i) for all large /i and A is such that the equilibration 
parameter is never small unless /i = avA, for some a at least order one. The analysis 
of asymptotics of restricted partitions in appendix ID] shows that for /i = ayK: 



r/^<cyAe5^('^)^, (4.31) 



where 




log(l - e-^) I - y^TT 

(4.32) 

and c is some order one constant. For A large enough so that £^ <^\, the equilibration 
parameter is small for all a < ajfi, where aj^ is the solution to g{ajii) = 0. 

Similarly, the dimension of the Hilbert space available to p^^ is 

A 

<^ = E^(^'^/^)' (4-33) 

where P{E, > p) is the number of partitions of E into parts at least p. Since these 
partitions have at most E/ p parts, then P{E^ > /i) < P{E^ E/p). It follows that p\^^ 
equihbrates if a > auv where auv is found by replacing a by A /a in (I4.32p and solving 
for g{auv) = 0. 

Higher-dimensional field theories 

In dimensions higher than one, it is easier to work with generating partition functions: 



kxj_,---,kx^ 
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For simplicity we assume that the field theory lives in a cubic box|^ Then, the free 
energy at large A is approximated by an integral: 



F{e-^)= --\ogZ{e-^) = n{d-l)p-^ / dk k'^-^ \og{l 
P Jo 

OO «QQ 

1 Jo 



g-/3fc^ 



p—fink 



= -r^'-^'^fid), (4.35) 

where f{d) = Q{d — l)r {d)({d + 1). The expectation value of energy and entropy are 
found in terms of (3: 

E = F + (5dpF = /3-(^+i) d f{d) 

S{(5) = (5^dpF = (5-'^ {d + 1) f{d). (4.36) 

The number of partitions at large E is well approximated by e'^'-^-': 

P{E) ^ e^^''''+'' (4.37) 

where A= (d+l) {d~'^ f{d)y^^'^'^^\ This is a well-known result for density of states in 
free d-dimensional theories. 

The partition function that generates P{E, s) is also easy to write down: 

\k\=s 
F, h h ^ ^ 

The free energy associated with this partition function is 



F =n{d-l)p-^y2 dk k'^-^ \og{l - e- 



(4.39) 



The expectation value of energy and the entropy are 

-^ = fi(rf-l)r(^+^)s('^+^)y^ ^ (4.40) 

S{P) = n{d - 1) (/3 ^' -^ - ^' dk k'-' log(l - e-^'^)) 



""^^Knowledge of the spectrum of Laplacian is sufficient to generalize this method to arbitrary compact 
spaces. 
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where we have defined t = [5k and ^ = /3s|lJ When £■ > 1 and ^i/li+i) > 1, ( I4.40P 
imphes that /3 is small. The entropy at small /3 is dominated by the term {s//3Y in 
f l4.4ip . Solving for entropy as a function of energy in this limit reproduces the familiar 
expression in f l4.37p for entropy of free theories. In this regime, since dg ~ o^a, the 
equilibration parameter is large. However, if ^^^d+i) = 0{1) then /3 is also order one, 
and P{E, s) grows as £''^/('^+^), but with a smaller coefficient. Similar to the situation in 
one spatial dimension, any subsystem of size s = aE^^^^'^^^ with a < ajR equilibrates, 
where aiR is the solution to g^am) = for 



g{a) = n{d-l)a'^ 2/3"^ 



^ dtf^ loKa-e-'^)\ A 



e*-l d 12 

„-.«,., „_i)/L_^|^. (,42) 

5 Atypical initial states and scrambling 

The typicality theorem implies that almost all initial states (typical states) in T-Lr 
equilibrate; however, in scrambling we are interested in far from equilibrium initial 
states which are highly atypical. For instance, consider an initial product state ip = 
4>M ® '^M for ^ small subsystem M. Then, from the fiuctuation theorem we know that 
for a fixed 4>m ^ 'Hm, and typical states ipM 



{\\Ps{t)-UJs{^MAM)\\)<\l , f . • (5.1) 

As a result, if d^j ^ d'g the density matrix of S is independent of ip^. As we saw 
previously, this condition is easily satisfied for small messageso However, the equi- 
librium density matrix can still depend on the message (pM- We say a Hamiltonian 
scrambles an ensemble of initial states, if all subsystems S smaller than half the size of 
the system us become independent of 0m- It is clear that this does not always occur. 
As was argued in [Ij , in systems where the eigenstates of the interacting Hamiltonian 
remain close to tensor products on M and M (as in some weakly coupled systems), us 
retains information about (pM- 

One might think that the type of arguments we have been using so far fails to 
provide insight about scrambling of small messages. The typicality arguments we have 
been using are based on the idea that a small subsystem S does not have enough 



^^For d — 1, after some algebra [4.401 reproduces the asymptotic behaviour we quote in appendix [Pj 
^® Replacing d^ with df^ in theorem |4. II implies typicality and universality for product states with 
typical tp eTifi . 
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Hilbert space (memory) available to keep track of all the information about the initial 
state. However, if we encode a small message (e.g. with one degree of freedom) in 
the initial state as in the case of a product state, it is likely that the message fits 
in subsystems with parametrically larger number of degrees of freedom. However, a 
careful look at theorem 14.11 shows that this is not necessarily true. Intuitively, if our 
Hamiltonian of interest is very efficient in generating entanglement, there might not be 
much room left in large subsystems to hold information. This is reflected in parameter 
5^ in (14. 6p . Consider the universality argument for initial state with ip]^}j fixed and a 
typical (pM e "Hm: 



\u;s{M)-{u;s)m\\)<\I^, (5.2) 

^^ = E^^It^I^) ^' [{trs{P''\k){k\P^))\ (5.3) 

where \k) are the eigenstates of the full interacting Hamiltonian, P*^ is the projector to 
"Hm- If the right hand side of (15. 2 p is small, the message M is scrambled in subsystems 
of size less than or equal to S. The parameter 6^ is a weighted average of regularized 
subsystems' purities in energy eigenstates. Note that if the subsystem purity (equiva- 
lently its second Renyi entropy) is bounded above (below) by an energy independent 
quantity a/dg, then 5^ dg < a. The logarithm on the right hand side of (15. 2 p provides 
a sufficient condition for scrambling: 

nM - log a >1, (5.4) 

where ua = \og{dA) denotes the number of degrees of freedom in subspace T-La- 

The Renyi entropy of subsystems in vacuum is a topic that has attracted a lot of 
attention recently. The ground state of systems with local interactions seem to have an 
entanglement entropy proportional to the area of the subsystem; see [TTj for a review. 
This goes with the name area law, and supports the simple picture that due to the 
locality of interactions, it is only the degrees of freedom near the boundary that matter 
in the calculation of entanglement entropy. The same is believed to be true for low- 
energy eigenstates and subsystems' second Renyi entropy. The sufficient condition for 
scrambling a message M in local systems with area law for all energy eigenstates in 
d > 1 dimension is: 

UM -ns-2c n^t^^''^ > 1, (5.5) 

and c is the constant appearing in the area law formula. There is some evidence that 
in discrete systems there is a high energy sector of the theory that has an extensive 
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entanglement entropy [T8| [T9] . In such cases, the coefficient c plays a crucial role in 
the size of the largest message that is scrambled. It is not clear whether such sectors 
survive in the continuum limit. 

With regards to the scrambling conjecture, the Hamiltonian of the matrix model 
or other quantum models of black holes include highly non-local interactions. For 
systems with non-local interactions it is not hard to imagine sectors of the spectrum 
that have extensive entanglement entropy. It seems to us that finding bounds on the 
scaling of entanglement entropy of energy eigenstates in matrix model or other models 
of quantum black holes is a promising approach to proving that they are scramblers. 



6 Conclusions and discussion 

In this work, we have generalized the recent information-theoretic approach to equili- 
bration to infinite-dimensional Hilbet spaces, including as examples field theories and 
matrix models. We generically find that small subsystems equilibrate, in the sense 
that their density matrices become almost independent of time and initial state. Our 
results suggest that at higher energies, larger subsystems can reach equilibrium. In the 
limit of infinite energy, we show that any subsystem smaller than half equilibrates. Our 
work emphasizes the importance of energy scales for equilibration and scrambling of 
large subsystems, both in field theory and matrix models. It is an interesting fact that 
even in Hermitian one-matrix model that has A^ degrees of freedom, the energy-scale 
relevant for equilibration of large subsystems is order A^^. This is in contrast with lat- 
tice field theories, where the relevant scale is the same order as the number of degrees 
of freedom N. 

For concrete estimation of the equilibration parameter we focused our attention on 
a quench problem from the high-energy sectors of free theories. However, the formalism 
is powerful enough to apply to more general ensembles of initial states. The challenge 
in each case is to find upper bounds on the dimensions of the Hilbert space available 
to subsystems. In our treatment of the problem, we find bounds on the support of the 
wave-function in the Hilbert space of free theory, and take advantage of our knowledge 
of the partition function at zero coupling. This is the reason why the dynamical 
constraint that appears is on the expectation value of the free Hamiltonian. Loosely 
speaking, the constraint we use, assumes an upper bound on the time average of particle 
number densities. We believe that this quantity is not a natural observable to consider 
in strongly coupled systems. It would be interesting to reformulate our approach in 
terms of energy expecation values in the interacting theory, and without any reference 
to the zero coupling Hilbert space, at least for some simple systems. 
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In spite of our powerful results in lattice field theories, we delibrately avoided ad- 
dressing the equilibration of density matrices corresponding to spatial regions in contin- 
uum field theories. The missing step in applying typicality to these cases is a faithfuo 
truncation of the reduced density matrix to a finite-dimensional Hilbert space. An- 
other subtlety appearing in the continuum limit is that one has to introduce a time 
scale the corresponding to turning on the interactiono However, this time-scale is 
clearly negligible compared to the time-scales required for the typicality argument to 
work. 

An interesting question to ask is whether one can generalize our discussion of equi- 
libration to multi-matrix models. In multi-matrix models, the notion of a subsystem, 
and how the true degrees of freedom interact is obscured by gauge transformation. Our 
philosophy is that for systems with gauge symmetry, the equilibration is more naturally 
formulated in terms of macro-observables discussed in appendix |B] [20] . 

Our focus in this work was on the equilibration of small and large subsystems 
by proving upper bounds on the time-averaged trace distance of states from their 
equilibrium value. However, in the context of information loss in black holes, one is 
also interested in deviations from equilibration. Finding a lower bound on the trace 
distance would be a natural way to address the question of how large is the smallest 
subsystem that does provably contain information about the full quantum state, but 
we leave this for future work. 
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Appendices 

A Levy's lemma 

Lemma A.l Given a function / : S''^ — ?■ M, with Lipschitz constant X = sup |V/|, at 
any random point x G S'^, 

Prob [fix) - (/),,e5. >e]< 2e-^('^+^)^'/^^ (A.l) 

where c = ttt- 

See [21] for a proof. 

A pure d- dimensional quantum state \ip) = J2k Cfcl^) can be thought of as a point on 
5*^"^"^ with the real and imaginary parts of c^ as coordinates. We apply Levy's lemma to 
functions f{ip) to discuss the probability of finding f{i/j) far from its ensemble average 

B Coarse-grained observables 

In a different approach to equilibration, Reimann [2] showed that, quite generically, the 
state of quantum systems becomes indistinguishable from its equilibrium with respect 
to coarse-grained observables. An observable is coarse-grained if it has finite precision, 
or in other words can have only a finite number of measurement outcomes. 

Following [3], we define the distinguishablility of two states p and a with respect 
to a set of coarse-grained observables Ai to be: 



D^(p,a) = i5^|tr(M,p)-tr(M,(T)|, (B.l) 



where M^ are all POVMs that describe the set Ai. The operational interpretation of 
distinguishablility is similar to the trace-norm: if Dj^ (p, a) < e the optimal probability 
of telling p and a apart using observables Ai is |(1 -|- e). 

The following theorem was shown in [3]. 

Theorem B.l Consider a closed quantum system evolving with a Hamiltonian which 
has non- degenerate energy gaps. The time-averaged distinguishablility of the state of 
the system from its time-average with respect to a set of observables M. satisfies: 

{DM{p{t)M)t<-A==,- (B.2) 
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where N{Ai) is the total number of outcomes in Ai. The denominator deff{u}) is 
exponentially large in the number of degrees of freedom. Therefore, for all realistic 
measurements we expect the state to be indistinguishable from its time-average. 



C Proof of theorem 14.11 : 



The proof we provide here for theorem 14.11 closely follows [T] . 

C.l Fluctuations: 

Proof of (i): Consider an initial state ipo = Ylk Cfc|^), where \k) is an energy eigenstate 
of the interacting Hamiltonian after the quench. Then, in the absence of degenerate 
energy gapscj, 

k,l 

= ^Cfcc;e-'(^^-^')*tr5 {P^\k){l\P^) (C.l) 

where we have used (^e~^^^''~^''')t = Ski- From the Cauchy-Schwarz inequality we know 
that the trace norm is bounded above by: 

(Wplit) - ^s\\)t < ^JdUtr[{p^s-^mt (C.2) 

The right-hand side of IC. 21 can be simplified further: 
{tr[{p^s-^int = E c,c„crc:(e-^(^-^'+^--^")*),tr[tr5(P'^|A:)(/|P^)tr5(P^|m)(n|P^)] 

= J2 \ck\'\ci\hr [trs (Pl^)(/|P^) tr 5 (P^|/)(A:|P^)] , 

= I]|cfc|'|Q|' Yl {sb\P^\k){l\P^\s'b){s'b'\P^\l){k\P^\sb') 

k^l s,s',b,b' 

= 5^|cfcnQ|' Yl {sb\P^\k){k\P^\sb'){s'b'\P^\l){l\P''\s'b) 

k^tl s,s',b,b' 

= 5^tr [tr5 (|c,pP^|A:)(A:|P^) tr 5 (|QpP^|/)(/|P^)] 

k^l 

= tr [{ulf] - J2 \ck\'tr [(tr sP^\k){k\P^r] 

k 

<tr[(u;|)2] (C.3) 



^^For a treatment of cases with degenerate gaps see [13] 
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where {\s)} and {\b)} are orthonormal bases that span Us and !-[§, and we have used 
the absence of degenerate energy gaps: (^e-^i^k-Ei+E^-E„)t^^ _ ^^_^^^^^_ xhe weak sub- 
additivity of purity [22] imphes tr [(cu^)^] < d^ti [{u^)'^]. Expanding in the orthonormal 
basis of free theory {|n)} it is clear that 

tr((u;'^)2) = J2 UJnmUJM{l\P''\n){in\P^\k) 

n,m,k,l 

El |2 

\^nm\ 

n,m,k,l e'HE<A 

< Yl \^nJ^ = tTitU^). (C.4) 

n,m,k,l 

Putting this back into ( 1C.2P we find: 



Mil 

deffiuj) 



Ip^sit) - u;l\\), < J -^^ = 4 (C.5) 
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C.2 Truncation Error 

Proof of {ii): Expand the state ?/'(t) in the orthonormal basis of free theory: 

V^(t) = X;Cn(t)cn.(t)|n)(m| (C.6) 

n 

with |n) denoting the eigenstate |ni, . . . , un)- The expectation value of the free Hamil- 
tonian at time t satisfies: 

oo oo 

>AtT{^-^^). (C.7) 

The matrix ip(t) = I , ^^^ J is positive semi-definite, and as a result ||A|P < 
II ^'^ II II V''^^ II ^ IIV^^^II [HI- From ()C.7p and the triangle inequality we find 



ll^(t) - ^p\t)\\ < 2\\A\\ + ll^^^ll < SyWn < 3^I^^^S 
Hence, the truncation error is bounded above by 

{\\{ps{t)-u:s)-{pl{t)-iol)\\), <(||(^(t)-^^)||),+ ||u;-a;^|| 
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where we have used the fact that partial trace as a quantum operation brings density 
operators closer, i.e. \\ps — csW ^ \\p — o'W- 

C.3 Typicality 

Proof of {in): We are interested in finding the probability for df.ff{uj) to be small. 
Let us define the following two maps: 

np] = Y.\k){k\p\k){k\ 
k 

m = Y.\k){k\p\k){ki (C.9) 

k 

where \k) are eigenstates of the interacting Hamiltonian, and \k) = |jp'^L||| are their 
projections to the initial ensemble. The first map F dephases states by killing off- 
diagonal elements in the interacting Hamiltonian basis, e.g. oo = ¥[iIjq\. The second 
map ¥[iIj] acts on states in T-Lt to dephase and projects them back to the initial ensemble 

An important observation is that the purity of ¥[iIjq] is larger than that of u, and 
therefore deffioj) > 1/tr (F['?/;o]^). Then, to prove {iii) it suffices to show that 

Prob ( tr (F[V^o]^) > jA < 26"^^^ (CIO) 

by applying Levy's lemma to the function 

/(^o) = ln(tr(F[^o]')). (C.ll) 

We need to find upper bounds on the average and the Lipschitz constant of /. We 
employ following lemmas: 

Lemma C.l For ip G I-Lr: 

( pA ^ pA\ /TT I (^\ 

(|^o)(^o| ® |^o)(^o|)^o = ^ ^^(^^ + 1) • (C-12) 

The swap operator § is defined by §|/, k) = \k,l), and we have used the trick tr (AB) = 
tr(§(A®5)). 



Lemma C.2 The Lipschitz constant of f{ijj) = In ( tr ( (F| ?/')(?/' |)^ j j is bounded above 
according to: 

A = |V/|<4rff. (C.13) 
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The proofs appear in Appendix B of [T]. With above lemmas in mind, it is not hard 
to come by an upper bound on (/)^o- 

{fU, <ln(tr(F[V^o]'))^o 

<lntr f§(F®F(|V^o)(V^o|)^o; 



dA{dA + 1) ^ J 



<ln 



<ln(A|, (C.14) 



Plugging (IC.13|) and flC.14p into Levy's lemma for /[V'o] gives: 



Prob Ttr (F[^o]') > ;^ ) ^ 26""^^, (C.15) 

which completes the proof of typicality. 

C.4 Universality 

Proof of (if): We apply Levy's lemma to the function gli^o] = \\uJs ~ i^s) -tpoW- "^^^ 
Lipschitz constant of g satisfies {Vgl < L As for the upper bound on {g)^^ we use the 
Cauchy-Schwarz inequality: 

< ^4{trM-{u;IU,y]h, (C.16) 
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The right hand side of IC.16l is bounded by: 



tr 



tr 



ss 



ss 



§ ( tr 55 
§ ( tr 55 



P^(F)®P^(F) (^o®^o)v., 






< tr 



ss 



< 



E 

Ik 



§ ( tr55A 

{kl\P^ ® P^\lk) 



p^m ® p^m 



p^m®p^m 



dA{dA + 1) 
(P^ ® P^)§) 



dl 



dl 



dl 



tr 55 [S (tr 55 [P^ ® P^|A;/)(A;/|P^ ® P^])] 



E ^^^''^^^''^'^^^ tr5 [tr5[P^|fc)(fc|Pltr5[P^|/)(/|Pl] 



Ik 



dl 



< E ^^'^^^l^^'^^'^^ tr5 [(tr5[P'^|fc)(fc|Pl)^ + (tr 5[Pl/)(/|Pl)^] 
Ik ^ 

d^ 



^Y.(k\p'^\f')^'s[{tTsP''\k){k\p^r] = ^ 



r/2 



(C.17) 



where we have defined P^{¥)[a] = P^F[a]P^ and used EH Therefore, {f)^, < 
^/d^5^^^d^. Plugging these into Levy's lemma for g{'4)Q) completes the proof. 



D Asymptotics of restricted partitions 



The large E asymptotic formula for the number of partitions of integer E is given by 

1 



P(P) 



4V3P 



r^2B/3 



(D.l) 



which is the first term that appears in the asymptotic expansion obtained by Hardy 
and Ramanujan in [23] . 

The asymptotics of the restricted partition function Q{E, s) defined by the number 
of partitions E into exactly s parts is studied [2l], and has the following three regimes: 



1. If s = 0(l) andP> 1: 



Q{E,s) 



E' 



s\{s-l)\ 



2. Ifs//E = 0(l) andE>l: 

a2 exp(f /;^-slog(l-e-")-f) 



Q{E,s) 



27rs2 



((e-i)/;#^) 



1/2 
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where a is found by solving 

E 1 r dtt 

/o 
3. If s/VE > 1, and E > s > 1: 



s^ a^ In e* — 1 



^/^ X 1 r^E Txs V6E _^ 
Q{E, s) ~ — ^— exp TTA/ -= e ^^ 



AVSE y \ 3 y/OE vr 

The number of partitions of E into exactly s parts is equal to the number of partitions 
of £■ — s into at most s parts: Q{E, s) = P[E — s, s). Since we are only interested in 
the limit E ^ s, the asymptotic behaviour of Q{E, s) and P{E, s) are the same. In 
section l¥^ we provide a physicist's derivation of these formulae. 

In number theory, the denumerant D{E]ai, . . . ,as) denotes the number of ways 
one can partition a positive integer E into integer parts ai, . . . , a^: E = ^^ JT-iCtj. The 
asymptotic behvaiour of the denumerant for large E and s = 0(1) is given by 

D{E; ai, . . . , a,) = ^ . (D.2) 

For a reference on asymptotic properties of denumerant see [25] . 



E Equilibration time 

In this appendix, we would like to elaborate on time-scales relevant for weak equilibra- 
tion. One way to approach this problem is to bound the trace distance ||p5(t) — ujs\\ 
averaged over a finite time interval [0,T]. This approach was discussed in [13], and a 
generalization of the fluctuation theorem to finite-time averages was obtained. Applied 
to our regularized Hilbert spaces we find: 



{\\pUt)-u^l\\)T<vl^l+'-^ (E.l) 

where (/)r = ^ /q dtf{t), d^ is the dimension of the total Hilbert space ip^ explores, 
and e is defined by the minimum spaces between energy gaps: 

e = mm{\G(a,b) - G(^c,d)\ '■ ("^ ^) 7^ {c,d)}. (E.2) 

a,b 

The minimum is over all energy gaps G(a,b) = Ea — Eb in the spectrum of the interacting 
Hamiltonian. We expect e to be at least as small as 0((i^^), which gives an equilibration 
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time-scale exponential in the number of degrees of freedom. This time scale is extremely 
long, and does not seem to provide any insights into the dynamics of equilibration. We 
could have anticipated this based on the fact that it takes a long time for ilj{t) to explore 
the exponentially large Hilbert space needed for the typicality argument to work. A 
stronger bound would require restriction to a specific set of Hamiltonians or further 
dynamical assumptions. 
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